Terrace- Width Distributions and Step-Step Repulsions on Vicinal Surfaces: 
Symmetries, Scaling, Simplifications, Subtleties, and Schrodinger 
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For more than three decades, measurement of terrace width distributions (TWDs) of vicinal crystal 
surfaces have been recognized as arguably the best way to determine the dimensionless strength 
A of the elastic repulsion between steps. For sufficiently strong repulsions, the TWD is expected 
to be Gaussian, with A varying inversely with the squared variance. However, there has been a 
controversy over the proportionality constant. From another perspective the TWD can be described 
as a continuous generalized Wigner distribution (CGWD) essentially no more complicated than a 
Gaussian but a much better approximation at the few calibration points where exact solutions exist. 
This paper combines concisely the experimentally most useful results from several earlier papers on 
this subject and describes some advancements that are in progress regarding numerical tests and in 
using Schrodinger-equation formalism to give greater understanding of the origin of the CGWD and 
to give hope of extensions to more general interaction potentials between steps. There are many 
implications for future experiments. 

PACS Number(s): 05.40.+j,61.16.Ch,68.35.Md,68.35.Bs 



I. INTRODUCTION 

Quantitative measurement of the widths t of terraces 
on vicinal surfaces became possible a decade ago. A 
principal motivation for examining the terrace width dis- 
tribution (TWD) is the recognition that it provides ar- 
guably the optimal way to assess the strength of the elas- 
tic (and/or dipolar) repulsion between steps, specifically 
the coefficient A of the elastic repulsion per length A/£ 2 . 
Here the elastic repulsion is taken to be perpendicular 
to the mean step direction. All standard analysis pro- 
cedures make a continuum approximation in the direc- 
tion along the steps, called y in "Maryland notation." 
(The perpendicular direction in the terrace plane, in the 
"downstairs" direction, is denoted x.) Hereafter, A ap- 
pears only in form of a dimensionless interaction strength 



A = Af3(k B T)- 



(1) 



where j3 is the step stiffness. 

Experimentally, a TWD is typically characterized by 
its variance a 2 and, at least when A is not small, has a 
shape that can be satisfactorily approximated by a Gaus- 
sian. The Gaussian form can be readily derived from a 
mean-field (Gruber-Mullins) QJ^] argument, which pro- 
duces an expression relating the variance to A. In recent 
years, theories from two new viewpoints have deduced 
different relations of A to the variance of the Gaussian. 
More recently, we have recognized that the TWD might 
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better be described using a simple expression arising from 
random matrix theory, called the "generalized Wigner 
surmise." As these results emerged, they have been pub- 
lished in several different articles The goal of the 
present paper is to collect succinctly the important re- 
sults, to provide a global view of progress on this prob- 
lem, to preview forthcoming results M , and to point out 
areas where further progress is needed. 

The following initial comments indicate our guiding 
philosophy: 1) The continuum approximation noted 
above is part of the step continuum approach to vicinal 
surfaces. In this perspective ||, the mesoscopic behavior 
of the step is characterized in terms of three parameters: 
the step stiffness $, the interaction strength A (or its 
equivalent), and a parameter representing the dominant 
kinetics (a kinetic coefficient or diffusion constant times 
carrier density). Hence, a knowledge of A is crucial to 
a proper description. 2) In this approximation, because 
step overhangs are physically forbidden, the set of step 
configurations in 2D space maps into the world lines de- 
scribing the evolution of non-crossing particles (spinlcss 
fermions or hard bosons) in ID space. This mapping is 
what leads to most of the progress in theoretical under- 
standing. 3) In experiments to date, investigators have 
measured the distribution of terrace widths £. This corre- 
lation function in essence is a many-particle correlation 
function, since one measures the probability of finding 
a pair particles separated by £ with none between them. 
(It is much easier for theorists to compute the probabil- 
ity of finding a pair particles separated by £, regardless 
of how many particles are between them; note that this 
two-particle correlation function should be equivalent to 
the many-particle one for step separations much smaller 
than the mean separation (£).) 

In experimental systems (cf. Table 2, below), A is typ- 
ically between and 15 0^,0] • (Whereas occasional val- 
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ues up to nearly 4000 have been reported S for A, our 
belief is that values above about 20-25 are indicative of 
anomalous behavior of some sort.) Exact theoretical in- 
formation is available only for ^4 = and A — 2 E^O , 
as well as in the limit A — > oo [p~2|— p~4[ . Hence, to as- 
sess the merits of various approaches for general A, we 
have generated well-characterized distributions numeri- 
cally. We have then compared each of the theoretical 
predictions with these calibration standards. 

In Sec. H we collect and synthesize the main results 
first for the traditional Gaussian analysis of TWDs and 
the competing ways of interpreting their variance in 
terms of A, then for the generalized Wigner distribu- 
tion arising from the theory of fluctuating systems. Sec. 
Ill recounts concisely several highlights of previous ex- 
plorations of these ideas. These include useful results on 
fitting procedures, an estimate of when discreteness be- 
comes important, and a procedure to gauge how many in- 
dependent measurements are contained in an image. Sec. 
[iv| gives a brief summary of findings in applications to 
experimental data, with an emphasis on trends. In Sec. 
^ we present previews of unpublished results concern- 
ing new directions in understanding TWDs with greater 
insight and in more complicated situations. Finally we 
offer brief conclusions and comments on connections with 
other active subjects in condensed matter physics. 



II. KEY RESULTS 

A. Gaussian Approximations to TWDs 

It is convenient and natural to divide £ by its aver- 
age value, thus constructing the dimensionless parameter 
s = £/{£). Then the TWD, P(s), is not just normalized 
but has unit mean. The Gaussian approximation to the 
TWD is then written: 



P(s) « P G (s) 
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(2) 



Gaussians are typically chosen, not just for their sim- 
plicity, but because their use can be justified readily for 
strong elastic repulsion between steps. In this limit the 
motion of each step tends to be confined near its mean 
position, a Gruber-Mullins (GM) argument (in which a 
single step is treated as active and its two neighbors are 
fixed at twice {£)) shows that 



48 in Kcm(NN) is replaced by 87r 4 /15 w 52, decreasing 
the variance by a scant 3 + %; i.e., KoM(all) ~ 0.139. 

The Grenoble group (ljjl^l pointed out recently that 
the variance in Eq. ([}]) using K^ai underestimates (for 
given A) the true variance. Their arguments are based 
on two ideas. First, the contribution of the entropic 
repulsion decreases with increasing energetic repulsion; 
physically, large energetic repulsions diminish the chance 
of neighboring steps approaching each other, where the 
non-crossing condition underlying the entropic repulsion 
becomes significant. Thus, for very large A the entropy of 
interaction becomes negligible, so that the only entropy 
is that of the individual steps. Secondly, if both steps 
bounding a terrace fluctuate independently, then the 
variance of the TWD should be the sum of the variances 
of the fluctuations of each step, i.e. twice the variance 
obtained in the Gruber-Mullins picture (in which there 
is a single "active" step between a pair of straight/rigid 
neighboring steps). This factor is reduced modestly by 
corrections due to the [anti] correlations |l7j of neighbor- 
ing steps. As a result, in this perspective the factor of 
48 in Kqm(nn) should decrease to 14.80, increasing the 
variance for a particular A by a factor of 1.801. 

Including entropic repulsions in an average way 
(mnemonically denoted X — EA, the two highlighted let- 
ters) H rather than discarding them extends to smaller A 
the range of viability of this (modified) asymptotic limit. 
Explicitly, A is replaced in Eq. (||) by an effective inter- 
action strength A e g obtained from the cubic term of the 
expansion of the projected free-energy of a vicinal surface 
as a function of misorientation slope. |l8| The resulting 
enhancement is 



A 



-2*=-L[ V4.1-I 
A 4A 



l+i~5 



A- 



(4) 



Explicitly, Eq. (||) becomes a 2 ~ KeaA~^ 2 , with val- 
ues for Kea given in Table 1. In this case, K ea(n n) is 
nearly 10% larger than -f^EA(a;z) -Q Since this approach 
represents the limit of minimally important entropic in- 
teractions, presumably this large ratio is an upper bound, 
approached for large A, while the smaller (3 + %) ratio of 
the GM case is more appropriate for weaker A. (In the 
free-fermion limit (A = 0), the NN and "all" cases obvi- 
ously must be the same!) 

The preceding approaches make a continuum approx- 
imation along the "time-like" y-direction but maintain 
discrete steps. By making a continuum approximation in 



r 2 = KxA- 1 ' 2 , 



(3) 



where the subscript X anticipates that there will be differ- 
ent proportionality constants in different approximation 
schemes, indicated by X. For the Gruber-Mullins case, 
with interactions only between nearest-neighbor steps, 



K, 



GM(NN) 



1/V48 w 0.144. For the Gruber-Mullins 
case, if all steps are allowed to interact with A/£ 2 , then 



*A value equivalent to K E A( a ii) — 0.247 was found explicitly 
in a calculation using the harmonic, lattice approximation of 
the Calogero-Sutherland model [n4|, as well as implicitly in 
earlier studies |l2Ul9| , po| | , and seems to be the exact asymp- 
totic coefficient fllll. 
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the x-direction as well and invoking correlation functions 
from roughening theory (so denoted X = R) , the Saclay 
group 

fl7fl2l[@ arrived at a result of the form of Eq. (§), 
again with A c g replacing A, in which Kr = 2/tt 2 sa 0.203. 

Since the various Gaussian approaches make different 
fundamental approximations, the detailed relationships 
between the width of the Gaussian and A differ notably. 
Even when a TWD can be well fit by a Gaussian, the 
estimation of A can be ambiguous. 



B. Symmetry and Wigner Approximation to TWDs: 
Continuous Generalized Wigner Distribution 

In considering high- lying energy levels in nuclei, 
Wigner long ago proposed that fluctuations in their spac- 
ings in energy should exhibit certain universal features 
depending only on the symmetry — orthogonal, unitary, 
or symplectic — of the couplings. This work, embedded 
in random-matrix theory ]10| , pTj| , has had profound and 
widespread implications for characterizing a wide range 
of fluctuation phenomena Jl^Jl^]. Since TWDs are an 
example of equilibrium fluctuations |^-^|, this body of 
knowledge should be applicable to them. The explicit 
connection is based on the description of steps using 
the (Calogero 23 1-) Sutherland |^4],^5| model of spinless 
fermions on a large ring (essentially ID with periodic 
boundary conditions) interacting with a repulsion decay- 
ing as the inverse square of separation. Remarkably, the 
distribution of interparticle spacings along the ring (i.e. 
the TWD) is equivalent to the distribution of the above- 
mentioned energy spacings, which can be solved exactly 
by random-matrix methods for the three symmetries fj] 
According to the so-called Wigner surmise, these three 
exact solutions for the distribution of fluctuations can be 
approximated by M 



The variance of P e (s) is just 



a a s' 



' exp (- 



(5) 



The three symmetries correspond to the values g = 1, 2, 
or 4, respectively. The constants b Q (associated with unit 
mean of P(s)) and a e (deriving from normalization) are 

2 



and 



2b 



(e+i)/2 



r(£±i) 



(6) 



This surmise was used to describe the spacings not of 
particles in real space but rather of energy levels, first in 
nuclei, later in chaotic systems fl26|j. 



^By construction, the wavefunction ^ of the Sutherland 
model is the product of the q/2 power of the differences of 
all pairs of particle positions, so that | >3> 2 1 was recognized to 
be identical to the joint probability density function for the 
eigenvalues of random matrices from a Gaussian |m| or a cir- 
cular P5l ensemble. 
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2b r 



(7) 



The approximations prove to be outstanding, accurate 
to better than ±0.004 for the latter two cases (cf. Fig. 1 
and Ref. [^6), Fig. 4.2a). From the mapping of the step 
problem onto the Sutherland Hamiltonian comes the 
relation 



A = q{q-2)/A, 



(8) 



(By inverting Eq. (||) to obtain g as a function of A, we 
can make the important and useful identification of g as 
2^/ieff, as given in Eq. (|)!) 

For the three special values of g, Eq. || accounts for 
the cases A = -1/4, 0, or 2, respectively. The value 
corresponds to steps interacting only via the entropic re- 
pulsion, whereas the negative value corresponds to an at- 
traction, which cannot be produced by the generic elastic 
interaction between steps (except perhaps in abnormal 
cases in which there are strong in-plane dipoles at the 
step edges |l^^]). The third case, A — 2, corresponds 
to a rather moderate repulsion. As documented in Ta- 
ble 1, the variance of Wigner's P e (s) is nearly the same 
as the exact value. The Saclay and the GM estimates 
are a few percent too low, while the modified Grenoble 
estimate is much too high. 

The crucial question is what to do for more general 
values of A. We simply use Eq. (^) for arbitrary value of 
g > 2, with g related to A by Eq. (||). For brevity, we 
refer hereafter to this distribution, for general g, as the 
CGWD (continuum generalized Wigner distribution). In 



P(s) 




FIG. 1. P(s) vs. s = £/{£} for the [sixth approximant Q 
to the] exact "free-fermion" , A—Q result (solid curve), the 
Gruber-MuUins approximation sin 2 (its/ 2) (long-short dashed 
curve), and the q=2 Wigner surmise result (dotted curve), 
barely distinguishable from the exact result). Offset upward 
by 0.4 for clarity, a similar plot of an approximant of the exact 
result for A— 2 J37J, the Gruber-Mullins Gaussian approxima- 
tion (24/7r 2 ) 1/4 exp(- v / 24(s - l) 2 )(long-short dashed curve), 
and the g=4 Wigner surmise result (dotted curve) . 
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contrast to the three special cases, there are no symme- 
try arguments to justify the CGWD form. We offer the 
following arguments in its support, although ultimately 
one must rely on numerical checks. 

1) It seems plausible that P g (s) is a decent approxima- 
tion of the TWD for values of g between 2 and 4 since 
the range in parameter space is small. In any case, the 
arguments supporting the approaches leading to any of 
the Gaussian approximations fail in this regime. 

2) Extrapolation of the CGWD to values of g greater — 
possibly much greater — than 4 is of more concern. For 
very large A, the argument underlying the Grenoble 
viewpoint becomes compelling. In this limit, the lead- 
ing term in the expansion of in Eq. (R) implies that 
Kw = 1/4 in Eq. (||), with A e g replacing A. Thus, as 
listed in Table 1, the CGWD variance approaches the 
[modified] Grenoble estimate nicely, while the Saclay es- 
timate is notably too small. Since the CGWD does well 
in the limit of very large g as well as at g=4, it is a promis- 
ing candidate for an interpolation method between these 
values. 

3) As a function of s, the CGWD not only has the Gaus- 
sian behavior expected (based on analogies with random 
walkers) at large step separations but also reproduces the 
exact power of s for s<l: In this limit, the many-step 
correlation function becomes identical to the pair corre- 
lation function, due to the vanishing probability of any 
other step lying between the pair of steps separated by 
s. Several workers have shown that in this limit, the pair 
correlation function is proportional to s B , with a pref- 
actor similar (within at least a few percent for physical 
values of g) to a e |l2]|^,|2§f(| . 

4) We can derive the CGWD from a Schrodinger- 
equation approach j(| , as discussed below in Section V B . 
This approach has the further benefit of allowing one 
to consider more general potentials than the asymptotic 
form of the elastic repulsion. 



C. Preliminary Numerical Results 

To test numerically the accuracy of Eq. we apply 
standard Monte Carlo methods to the most elementary 
model that contains the necessary physics, the terrace- 
step-kink (TSK) model. In the TSK model the only ther- 
mal excitations are kinks of energy e along the steps. 
The stiffness (3tsk of an isolated step — needed to ex- 
tract A from A— is simply 2k B T(a\\/a 2 L ) sinh 2 (e/2fc B T) 
[pip . Here a|| is the unit spacing along a step edge (y), 
and aj_ is the x component of a kink. This model is 
obviously discrete in the y as well as the x directions 
P, p7|j2^j3l[ | . For simplicity we consider a vicinal simple 
cubic lattice with unit lattice constant: a\\ — a± = 1. 
Periodic boundary conditions are imposed in both di- 
rections. To minimize finite-size effects, the length of 



the system in the y-direction, L y , should be substan- 
tially larger than the characteristic distance y co \i along y 
between close approaches of adjacent steps: (£) 2 /3/kBT 
pl| - The choice of the mean spacing between steps re- 
quires particular care. We shall show below that if (£) is 4 
or smaller, finite size effects may contaminate the results 
extracted from the CGWD (since it is based on a contin- 
uum approximation). On the other hand the minimum 
acceptable value of L y increases like (£) 2 . Furthermore, 
too low a temperature results in slow dynamics and a 
high stiffness, making demands on L y , while too high a 
temperature leads excessive step wandering and break- 
down of the approximations underlying the viability of 
the TSK model. We are preparing a careful discussion of 
these considerations @, which includes transfer- matrix 
calculations in addition to Monte Carlo simulations. 

In Fig. 2 we provide some preliminary results for the 
case (£) = 6 at fcsT/e = 0.5, with L y — 200 and the num- 
ber of steps N = 10. Data from runs using the standard 
Metropolis algorithm |52| were taken, after equilibration, 
over 10 5 Monte Carlo steps per site. In addition we use 
the "refusal-free" n-fold way p3p^|, especially for large 
A (or at low T). There it is much more efficient than 
the Metropolis algorithm, which requires many attempts 
before making a change. The elastic repulsion is here 
considered only between neighboring steps, a common 
simplification in Monte Carlo 0|2l]], with the accordant 
modest underestimate of a 2 noted in Subsection A. (In 
Ref. we will also extend the inverse-square repulsions 
to further neighbors.) Our algorithm includes "corner 
exclusion" in addition to standard edge exclusion, based 
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FIG. 2. Plot of the variance a as a function of A on a 
logarithmic scale, plotted for the CGWD ["Wigner distribu- 
tion" ] (solid curve) , the modified Grenoble (dotted curve) and 
Saclay (long dashed curve) Gaussian distributions, and the 
Gruber-Mullins Gaussian approximation (long-short dashed 
curve). The CGWD curve passes essentially directly over the 
exact value of the variance at A=2. Monte Carlo data gener- 
ated using the Metropolis algorithm are depicted by Q's; data 
produced with the n-fold way algorithm are shown as »'s. 
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on some evidence that it provides the better discrete ana- 
logue of the continuum model; the consequent modest up- 
ward shift of cr 2 is in the opposite direction of that due to 
restriction to nearest-neighbor step-step repulsions. (See 
Ref. § for more details.) 

Along with the numerical results, the various predic- 
tions of the variance are plotted as functions of the phys- 
ical variable A. A logarithmic scale is used for the hor- 
izontal axis so as not to give undue visual emphasis to 
larger values of A nor to blur the region of rapid vari- 
ation for small A, for which an exact calibration point 
exists. The physical values of A range from near up to 
the mid teens. A few larger values have been reported, 
but there are suspicions that more than simple elastic 
repulsions are involved. There are relatively few reports 
of small but non-zero values of A. We suspect that one 
reason is that any of the Gaussian approximations mani- 
festly fail in this regime because the distribution becomes 
strongly skewed.^] Before the recognition of the utility 
of the Wigner distribution, one could not deal quantita- 
tively with small A (p5[ . 



III. USEFUL RESULTS FOR INTERPRETING 
EXPERIMENTS 

A. Extracting A from Variance 

If one accepts the CGWD as the optimal way to ana- 
lyze TWDs, then Eq. ([?]) shows how to estimate the vari- 
ance from A. However, experimentalists usually seek the 
reverse. An excellent estimate jj| of Aw from the vari- 
ance can be derived by expanding a w as given in Eq. (Q) 
in powers of This scries can then be reverted to give 
g as a function of a 2 ||. Then using Eq. (||) gives the 
estimate 



A 



w 



1 

16 



( CT 2)-2 _ 7 (a 2 )- 1 + 



27 
T 



:!r> 



(9) 



with all four terms needed to provide a good approxi- 
mation over the full physical range of A. The Gaussian 
methods described earlier essentially use just the first 
term of this expression and adjust the prefactor. When 
A is not weak (see Ref. || for explicit guidelines.), those 
who for some reason prefer not to use Eq. (||) to gauge g 
(and thence A) can extract the variance from a Gaussian 
fit and then applying Eq. (JsJ) is a reasonable procedure. 

When dealing with tabulations of data analyzed in the 
traditional way fq], i.e. using the inverse of Eq. (p|) with 



'"However, the idea of estimating A using this skewness || 
did not prove to be fruitful when confronting experimental 
data 101. 



X = GM(NN), it is useful to recast Eq. @ in a form 
that indicates the factor by which the estimate Aw based 
on CGWD exceeds the traditional estimate Agm(NN) 
(denoted Aq for brevity): 



A W /A G [= A w /Ac 



3 - 21a 2 



+ f- 6 • do) 



As noted parenthetically, Eq. (|J) implies that the ratio of 
the physical interaction strengths is the same as that of 
the dimensionless strengths. Since A cx A we can use this 
relation in Table 2 to update tabulated A's in Ref. . 



B. Gaussian Fits of the Generalized Wigner 
Distribution 



Since TWDs for strong repulsions are well described 
by Gaussians, one expects — and finds — that the CGWD 
should be well approximated by a Gaussian in this limit. 
In Ref. H a quantitative assessment is given of how 
closely the two distributions correspond as a function 
of g. At the calibration point (for which an exact so- 
lution exists) for repulsive interactions (g = 4), the rel- 
ative difference of the standard deviation of a Gaussian 
fitted to P e (s) from the actual standard deviation of this 
CGWD (viz. the square root of the second moment of 
P B {s) about its mean of unity) is around 1%, and de- 
creases monotonically with increasing g. For this range 
(g > 4) differences between estimates of A obtained from 
CGWD and the various Gaussian fit methods are pre- 
dominantly due to different philosophies of extracting A 
from a rather than from differences in the fitting meth- 
ods. 

In contrast to the Gaussian approximations, the peak 
of the CGWD must perforce (due to unit mean) lie below 
one. Specifically, for A = and 2, the maximum of P g (s) 
occurs at s = 0.886 and 0.940, respectively, while the 
limiting value for strong repulsions is 1 — 0.125A/A c ff [0- 
Formulas have been derived (S) indicating the errors in 
fitting g due to errors in the first or zeroth moment of 
the distribution. 



C. Wigner Distribution as a 2-Parameter Fit 

In applications to experimental TWDs, the CGWDs 
giving the best fits sometimes have first moments that 
differ somewhat from the first moments of the data, es- 
pecially in cases termed "poor data" which exhibit 
a small "hump" at large values of s, beyond the peak 



near unity [see Section IV below]. Moreover, it can be 



desirable to determine the scaling length (the "effective 
mean," which equals the first moment for ideal CGWDs) 
and the variance in a single fitting procedure rather than 
to predetermine this length from the first moment. This 
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"refined" scaling implies that the argument of P g should 
be £/£, where £ denotes the characteristic length deter- 
mined along with g in a two-parameter least-squares fit 
of the data to a CGWD. Since s is still determined from 
the raw data as £/(£), the refined scaling translates into 
replacing s by s(£}/£ in the argument of the distribution. 
If the integration variable s were similarly replaced, then 
the refined scaling would amount to a redefinition of a 
dummy variable, and normalization would still be real- 
ized. Since the independent variable is kept as s, we make 
the replacement: 



p e (s) - {{i)/i)P e W)/i) 



((£)/£)P g (£/£) (11) 



In the specific applications to data in subsections A and B 
of the next section, (£}/£ tends to be greater than unity, 
typically by several percent, but it is unclear whether 
this is true for semiconductors or other metals. In our 
companion Monte Carlo simulations JtJ, where we have 
greater control of purity and uniformity than in experi- 
ments, the optimal £ is essentially identical to (£): there 
is no need for the added flexibility of the two-parameter 
fit. 



D. Effects of Lattice Discreteness 

For actual crystals as well as for the TSK model used 
in numerical simulations, the variable s cannot assume 
a continuum of values as implicitly assumed in writing 
Eqs. (H|j); the only possible values of £ are integer mul- 
tiples of a±. If this restriction is placed on the values 
of £ used to generate the arguments s of Eq. (||), then 
we have constructed a discrete generalized Wigner dis- 
tribution (DGWD)J^] We use the same value of b e as in 
Eq. even though it is no longer guaranteed to pro- 
duce unit mean (or the same variance) as it does for the 
CGWD. Since these vicinals are technically rough, there 
is no need for (£) to be an integer multiple of a± (or oth- 
erwise in registry with the terrace plane), though it is 
common to make this choice in simulations. 

Scaling of discrete TWDs for the free-fermion case 
(A = 0) was demonstrated nearly a decade ago In- 
spired thereby, we [|| explored the effects of discreteness, 
first choosing values of (£) and g to specify a DGWD, 
then numerically performing two-parameter fits using 
CGWD formulae [Eqs. (|), (§), (0)] to produce esti- 
mates of g c and (via Eq. (@)) A c . 



^Discreteness also introduces the possibility of a roughening 
transition from a vicinal to a high-index- facet surface p2| , p6| . 
This interesting phenomenon does not seem to play a role in 
the physical systems under study here. 



Among many minor observations, two major themes 
stand out: First, (£) > 4, A c provides a reasonable es- 
timate of A over the range of physically reasonable di- 
mensionless repulsions. Furthermore, at fixed values of 
A the error in A c diminishes as (£} increases. Second, as 
the TWD becomes narrower (i.e. for sufficiently large A 
or q), A c becomes a questionable estimate for A; study 
of the cases (£)/a± = 2-6 suggests that this breakdown 
occurs for g near ((£}/a±) 2 . This threshold corresponds 
to the squared interstep spacing being comparable to the 



variance. 



For very large A, seemingly just above the range of 
greatest physical significance, there are more general in- 
dications of the breaking down of the continuum approx- 
imation. E.g., for A in the upper teens, there begin to 
be ambiguities in the application of Eq. (|§|) Q| . With pe- 
riodic boundary conditions one can still get elementary 
excitations that are extended along a step (i.e. along y), 
so long as they are "in phase" in x, but with more re- 
alistic conditions (with various sorts of defects hindering 
the fluctuations of occasional steps) , the elementary ex- 
citation becomes individual "teeth" (kink-antikink pairs 
separated by one spacing along y) |3§| ] . Then the idea of 
step stiffness also breaks down, and with it the concept 
of 1 (see Eq. (@)). 

The main implication is that analyses of highly mis- 
oriented vicinal surfaces with CGWD should be viewed 
with caution. E.g. the (1,1,7) for close-packed steps on 
surfaces vicinal to {1 0} planes of fee crystals corre- 
sponds to (£) = 3. For {1 1 1} fee surfaces, the corre- 
sponding Miller indices are (5 3 3) for A steps ({1 0} 
microfacets) and (2 2 1) for B steps ({1 1 1} microfacets) 

The obstacles posed by discreteness are not vagaries of 
Wigner distributions. High misorientation causes simi- 
lar problems when the mean and variance of discretized 
Gaussian TWDs are analyzed as though they were con- 
tinuous Gaussian functions. (See Ref. g for details.) 



E. Estimate of Number of Independent 
Measurements 

In order to estimate uncertainties in the determination 
of the TWD and, ultimately, A, it is important to have 
a realistic value of the number of independent measure- 
ments, a number generally much smaller than the total 
number of measurements. To make a rough estimate, one 
can compute the correlation function |40| of the terrace 
widths £ n (y) between steps n and n+1: 



C n (y) 



(N-n)(L y -y) 



{if 



(P) {If 



(12) 



is calculated, where N is the number of terraces in the 
image. The correlation function along the steps decays 



G 



exponentially as Co(y) ~ exp(— y/£ y ), where £ y is propor- 
tional to j/coii (cf. Eqs. (5), (12), and (26) of Ref. @), 
but can be measured directly. The correlation function 
between steps is more complicated. As noted in Table 1, 
Ci(0) is negative JhJ; |C„(0)| tends to decrease rapidly 
with increasing n. Setting c as a small cutoff (c = 0.1 is 
recommended Q), we determine y c , the smallest value 
of y for which |Co(y)| < c when y > y c , and n c , the 
smallest n so that |C„(0)| < c for all n > n c . Then the 
number of "independent" terrace widths will be approx- 
imately (L y /y c )(N/n c ) rather than L y N , as might be 
naively guessed. A rough test calculation Q shows that 
the reduction factor can be nearly two orders of magni- 
tude, emphasizing the need for using several STM images 
to obtain decent statistics. 

IV. APPLICATIONS TO EXPERIMENTAL TWDS 

In two papers j|,D , we made extensive applications of 
the ideas presented above to M. Giesen's voluminous data 
on vicinal Cu {100} and {111} surfaces, each at three 
different misorientations, and these six cases at various 
temperatures. In all, around 30 different cases were con- 
sidered. In addition, our ideas were tested successfully 
on data for vicinal Pt(110), which has a small A and so 
is not amenable to the Gaussian approaches used hereto- 
fore. The purpose of this section is to summarize the 
tabulations and discussions in those papers. 

A. Copper: Moderately Strong Repulsions 

The Cu TWDs can be sorted into three groups based 
on a visual assessment of their quality A "good" 

TWD changes height essentially monotonically below the 
peak and again above it; there is minimal scatter in the 
data points. An "OK" TWD has more scatter, with small 
dips and peaks introduced by variations (within the lim- 
its of the general margin of error) of single data points. 
A "poor" TWD has a double-peak or hump at large s; 
correspondingly, the position of the (main) peak occurs 
noticeably below s — 1, even when the peak is fairly nar- 
row and the skewness minimal. The judgment that this 
data is "poor" is based both on the intuition of the exper- 
imenter and on the following argument: A second peak at 
large s would be characteristic of the onset of faceting; 
however, "poor" data tends to occur at high tempera- 
tures, whereas faceting should be more important at low 
temperatures. 

The data fits exhibited several general trends. In al- 
most all cases, the value of £ derived from the two- 
parameter fit to a CGWD is smaller than (£) given by 
the mean of the TWD (and the opposite shift in the ex- 
ceptional cases is very small); likewise, the directly mea- 
sured values of the variance are almost always larger than 



the values obtained by any of the three fitted curves (cf. 
Sec. 7 of Ref. ||). The value of g is higher for the two- 
parameter CGWD fit than for the single-parameter ver- 
sion, and the associated value of a 2 typically closer to 
that deduced from the Gaussian fit. For "good" data, 
£/(£) differs from unity by a few percent, and the change 
in g and cr 2 is negligible. For "poor" data, 1/(1) is at 
least twice as far below unity, and the two-parameter- 
fit curve is narrower than the single-parameter-fit curve. 
The tails or humps in the experimental TWDs seem to 
be responsible for the systematic discrepancies in the fits, 
especially the smaller mean and smaller variance of the 
fits relative to the direct measurements. 

A remarkable consistency check was obtained for Cu 
(1 1 13) For a dozen values of temperature, (kBT) 2 A 
was plotted against T. Since A is expected to be rel- 
atively insensitive to thermal change, Eq. (|l|) predicts 
that the plotted curve should decrease like the stiffness. 
To within error bars, such behavior is found, where the 
stiffness is computed using an independently determined 
kink energy. 

B. Platinum: Weak Repulsions 

On vicinal Pt(110) at room temperature, the terraces 
are (1x2) reconstructed, and the steps correspond to 3- 
unit "(1x3)" segments. Recent measurements show that 
the interaction between their steps is small p5|, render- 
ing Gaussian approximations invalid. Fits to the CGWD 
yield g = 2.06 (A = 0.0309) or, when done in the two- 
parameter way, g = 2.24 (A = 0.134) ||; in the latter 
case, the optimal £/(£) is 91% and the fit is notably bet- 
ter. The presence of a high-s bulge indicates this feature 
is not peculiar to the vicinal Cu systems of GE. 

C. Other Systems 

Additionally, in Table 2 we list the variances mea- 
sured for several different experimental systems, along 
with the value of A deduced from the CGWD distribu- 
tion via Eq. (||). The primary goal of this table is to 
display general trends in physical systems rather than 
to provide a comprehensive account of experiments to 
date. As asserted earlier, values of A are generally be- 
low the mid-teens. On the other hand, A ranges over 
orders of magnitude. If one accepts that the CGWD pro- 
vides a good accounting in general for A as a function of 
the measured variance, then the column labeled Aw / Aq 
shows that for most systems, the underestimate by using 
the Gruber-Mullins approximation is roughly half that of 
the asymptotic limit. 
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V. NEW DIRECTIONS 

A. Multistep Distributions 

Experiments to date have focused exclusively on the 
TWDs, ignoring the possibility of extracting the distri- 
butions of the distances between pairs of steps having 
n steps, n — 1,2, or more, between them. This supple- 
mentary data could provide a valuable consistency check. 
For the three special cases p=l,2,4, these distributions 
have recently been investigated theoretically in a differ- 
ent context @. If in Eq. (|) we make the redefinition 
s = (£i + ... + £ n+ i)/(£) {I being the terrace width), 
then this CGWD expression gives a good approximation 
of the multistep distribution, provided that the power- 
law exponent g is replaced by 



(n+l)(n + 2) 
Qn = n H g, 



(13) 



where g (or, equivalently, g ) is the exponent for the 
[single] TWD. The new constant b Bn is determined by 
the condition that the first moment of P g (n, s) is n + 1; 
besides replacing g by g n in the T-function arguments in 
Eq. |6[ a factor of n+1 must be included in the denomi- 
nator. The normalization constant a 8rl can be obtained 
simply by using b Qn and g n in the expression for a e in 
Eq. (|(|). As for TWDs, these results can be taken to 
apply to general values of g. 

Preliminary checks using Monte Carlo simulations of 
the TSK model, described above, find fine agreement 
with this multistep CGWD for the double-terrace-width 
(n=l) case, but just adequate agreement for the case 
n=2. Moreover, Table 1 shows that the variance of the 
sum over the widths of two adjacent terraces predicted 
by P e (l,s) does not display the spectacular agreement 
with exact results seen for simple TWDs, viz. n=0. We 
suspect that the agreement will further degrade as more 
widths arc included (higher rt's considered), due to weak- 
ness in the main assumption in the derivation of Eq. ( 13 ) : 
that the conditional probability density of occurrence of 
a step at a given distance from a fixed step, with n steps 
in between, can be expressed in terms of the (n+l)th 
power of the corresponding probability for this distance 
with no intermediate steps. 



B. CGWD and Beyond via Schrodinger Equations 

As presented, the CGWD is formally justified only for 
the three special values of g. Accordingly we have de- 
veloped arguments using Schrodinger equations to show 
that it can be expected to have the more general valid- 
ity assumed above ||. The formalism also allows treat- 
ment of more general potentials than the inverse-square 



term characterizing the long-range behavior of elastic in- 
teractions. Of particular physical importance are the 
higher-order terms that enter at smaller terrace widths 
and an oscillatory interaction mediated by electronic sur- 
face states. 

We begin by defining a wave function ipo(s) such that 
ip 2 (s) = Pg(s). Differentiating twice and using Eq. (||), 
we find 



rV'o(s) 



As- 



b e (g+l) + b 2 e s 2 



Ms) = o. 



(14) 



The term V(s) = As~ 2 is the dimensionless step interac- 
tion. The term U(s) = b 2 s 2 is a dimensionless projected 
free energy representing interactions with all the other 
steps not explicitly considered. Clearly Eq. ([u]) can be 
understood as a Schrodinger equation, with ipo(s) the 
(real) ground state wave function and with b B (g + 1) the 
associated eigenvalue. (To consider perturbations from 
pure inverse square interactions, one can generate all the 
eigenf unctions ip n {s) of Eq. (fu]), which can be expressed 
as special functions (see Ref. || for details).) In this 
framework, by substituting more general potentials for 
As~ 2 in Eq. (|l4|) and solving for the ground-state wave- 
function, we can contend analytically with more compli- 
cated potentials. Successful tests are described in Ref. [||. 

It is tantalyzing to invert the preceding approach to de- 
duce the underlying interaction potential from the exper- 
imental TWD. Since naive implementations prove to be 
dangerous, the recommended procedure || requires con- 
siderable computation: Initially, a parametric approxi- 
mant of the unknown potential should be constructed us- 
ing all available information; crude initial estimates must 
be made of the values of the parameters. If the tail of the 
experimental TWD is Gaussian, a monotonically decay- 
ing V(s) and a quadratic t/(s) are anticipated. (If the tail 
is exponential, U(s) is linear in s, and V(s) can be non- 
monotonic.) Then the choice of parameters is optimized 
by iteratively minimizing the least-squares difference of 
the experimental TWD and |^q(s)| 2 , where V'o(s) is the 
numerical solution of Eq. (14) or its equivalent. In a test 
of this procedure, for very large or very small s the de- 
rived potential V(s) was significantly different from the 
potential used to generate the "experimental" TWD, but 
they were quite close to each other over the range where 
the TWD is large, 0.5<s< 1.5. Hence, to improve esti- 
mates of the potential over a large range of I, one should 
fit TWDs measured for several different misorientations 
(and, if possible, for different temperatures). 

Applications to vicinal Cu surfaces is in progress |fi2| . 
Preliminary fits can account for secondary humps men- 
tioned in Section [V by invoking nonmonotonically de- 
caying interactions and exponentially decaying tails. 
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C. Oscillatory Interactions Mediated by Surface 
Electronic States 



The preceding subsection began with allusions to oscil- 
latory electronic indirect interactions between steps [Esl . 
When mediated by bulk electronic states, such interac- 
tions between atoms on surfaces decay rapidly with sep- 
aration, but if mediated by circularly-symmetric surface 
states, the envelope of the oscillatory interaction has the 
same inverse-square behavior as the monotonic elastic, 
dipolar, and entropic repulsions: 



£~ 2 cos(2k F e + <f>), 



(15) 



where fc^ is the wavevector of the surface state at the 
Fermi level and ^ is a phase shift associated with scat- 
tering from the pair of steps. 

Consistent behavior was seen in measurements of 
TWDs on vicinal Ag(llO), including the presence of a 
surface state in the appropriate place in the surface Bril- 
louin zone [EI] . However, the evidence for the influence of 
surface states was not compelling due to the large number 
of fitting parameters compared to the amount of experi- 
mental data. 

Convincing evidence of long-range, surface-state medi- 
ated interactions between Cu atoms on Cu{lll} has just 
appeared . It is tempting to invoke these interactions 
(which could even decay as £~ 3 ^ 2 due to the isotropy 
of the state as the source of the large- s humps in 
"poor" data on this surface. However, this idea does not 
explain similar "poor" data on Cu{100}, where the image 
states are far from the Fermi energy. 

The energy of the long-range interaction measured for 
Cu atoms on Cu{lll} is notably weak [EBj: the deep- 
est minimum corresponds to an attraction of 0.4 meV at 
27A. However, the consequent prefactor of the expres- 
sion m Eq. (H|) is 0.3 eV- A 2 , which is comparable to the 
values of A listed in Table 2 (if some length of order an 
atomic spacing is used to adjust the units). 

The transport properties of fractional metal overlayers 



have received close scrutiny in recent years 1 47 . Since the 
metallic surface states can be tuned in these systems, it 
is intriguing to speculate about engineering morphology 
using such step interactions. 



VI. CONCLUDING REMARKS 

The CGWD of Eq. (|J) is an excellent interpolation be- 
tween the established points at A=0 and A=2, and ap- 
proaches the correct limit for very large A. Qualitatively 



"For non-circular states, the value of the wave vector in 
Eq. ([[5]) is that along the Fermi "surface" at which the elec- 
tron velocity is in the x direction; see Ref. p3| for details. 



it certainly captures the global behavior of variance as a 
function of A, and numerical evidence suggests that it in- 
terpolates well in the regime of large A. While the shape 
of the TWD does approach a Gaussian in this regime of 
moderately strong A, the CGWD (via Eq. (0) provides 
arguably the best way to extract A from the variance of 
the TWD. Of the several ways to extract A from fits to 
a Gaussian, the Saclay (R) scheme is better for moder- 
ate A while the Grenoble (EA(all)) scheme is better for 
stronger A. 

The difficulty of accurately estimating the value A from 
the TWD, especially from its variance or width, is exac- 
erbated by the extreme sensitivity of a 2 to A: fractional 
errors in the deduced widths of TWDs are magnified by a 
factor of 4 in A. For many applications, determining the 
relative size of the step-step repulsions between different 
systems is more important than deducing their absolute 
sizes; in such comparisons, it is crucial that the analysis 
of A be done using the same approach for all systems. 
(Likewise, experimentalists should state clearly the raw 
[dimensionless] width, a — or the value of g in a fit to the 
CGWD.) 

Often the extracted value is rationalized by misap- 
plying the celebrated result of Marchenko and Parshin 
relating the step repulsions to surface stress [ Egl . 
That formula assumes an elastically isotropic substrate 
and asymptotically large separations. Usually one, often 
both, of these conditions do not apply, and there is no 
well-prescribed procedure to compute corrections. Fur- 
thermore, the in-plane component of the stress dipole 
is not measurable and is often neglected. Thus, estab- 
lishing quantitative connections between deduced A and 
surface stresses is [even] harder than extracting reliable 
quantitative estimates of A. 

Another worrisome assumption is that the step inter- 
actions are "instantaneous" in the 1+1 D perspective (i.e. 
occur only between points on steps at the same coordi- 
nate along the mean step direction y) becomes particu- 
larly questionable when steps are close together and have 
large wandering fluctuations with short wavelengths. In 
the latter situation, the description of single-step fluctu- 
ations in terms of stiffness may also break down. 

In systems in which surface states near the Fermi en- 
ergy play an active role, there should be notable effects 
on the TWD and the consequent surface morphology. 
Multistep correlations have received little attention, even 
though data is readily available in experiments measur- 
ing TWDs. Moreover, it is almost as easy to tabulate 
the step-step pair correlation function as the TWD, but 
easier to decipher theoretically. We have provided several 
hints and warnings, hopefully useful, for experimentalists 
studying spacings on vicinals. 

Most theoretical activity dealing with random fluctu- 
ations in complex systems has focused on the the three 
special values of g, and occasionally on interpolations 
between them. This corresponds to weak A. The few 
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exceptions focus on pair correlations and are rather tech- 
nical [p~2| | . More explicit numerical investigations in this 
regime would be illuminating. Moreover, there remains 
the mystery of why the CGWD works so well when there 
is no fundamental symmetry argument to justify it. 

The Calogero-Sutherland model has been termed an 
ideal |l4j Luttinger liquid jl9 50 1 . Connections have been 
made to edge states in the quantum Hall effect |5l[] , non- 
linear waves in a stratified fluid ^2j, and a host of more 
abstract problems. Most of these systems exhibit cor- 
rections, making it difficult to make detailed connections 
with the theory . ft will be interesting to see whether 
similar problems involving corrections to the A/£ 2 much 
discussed above confound a similar effort for vicinal sur- 
faces. Furthermore, many of the interesting properties 
involve dynamic correlations, which for vicinal surfaces 
translates to correlations between displacements on [dif- 
ferent] steps at different values of y. In any case, however, 
these connections between the properties of vicinal sur- 
faces and other active fields add to the fascination of the 
subject. 
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Table 1: Tabulation of predictions of the variance of terrace- width distributions P(s) [where s is the terrace width 
normalized by its average value] based on exact results at the three soluble values of the dimensionless interaction 
strength, A, the corresponding Wigner-surmise expression, and several ways of interpreting a Gaussian fit. 



Property 



Case 



Abbrev. 
X 



1 + y/T+AA = 2^A eS 



Ref. 



Q = 2 

Non-interact 
1 = 



g = 4 
Exact rpl. 
A = 2 



Arbitrary g 
Repulsive 
A=(q-2)q/4 



g — > oo 
Extreme rpl. 
A^A cS = g 2 /A 



Symmetry assoc. w/ Sutherland TL 



unitary symplectic 



~5TT 



[SHO+phonons] 



32/tt 2 
4/tt 



(64/9tt) 3 
64/9tt 



In leftmost 
column 



2exp[(g + 3)/2] 
g/2 + 1/4 



Variance 
o 2 = /i2 
= /4-l 



Exact [all] 
Wigner surmise W[all] 



Gruber-Mullins 



Gaussian modified Grenoble 
alternatives " 

Saclay 



GM(all) 
GM(NN) 

EA(all) 
EA(NN) 
R 



f 



0.180 
0.1781 



0.105 
0.1045 



(e+i)/26 fl 



1 



0.495/p 
0.500/g 



Neighboring Exact ((si + S2 — 2)' z ) 
terraces Multistep CGWD 



0.1307 

0.1307 

0.247 

0.260 

0.203 



0.0981 
0.1021 
0.1185 
0.1300 
0.101 



0.139 /\M 
0.144 A/i 
0.247/Vieff 
0.260/V^cff 

0.203/Vieff 



0.278/g 
0.289/p 
0.495/g 
0.520/p 
0.405/p 



0.248 
0.257 



0.138 
0.145 



The "unmodified" Grenoble expression (T^,|T^] is obtained by substituting A for A\,g. The entries for the variance at 
A=2 are almost 50% larger. The bracketed "[all]" is a reminder that the Calogero-Sutherland model — and, hence, the 
Exact solutions and the Wigner surmises — involves all steps interacting. As A increases, the TWD becomes narrower, 
more symmetric, and more nearly Gaussian. Anticorrelations of neighboring terrace-width fluctuations increase. For 
the three exactly-solvable [non-trivial] cases, the Wigner surmise provides an excellent approximation, significantly 
better than any alternative. 

The covariance of the fluctuations of neighboring steps is the difference from unity of the tabulated correlation 
((sj + s 2 — 2) 2 ) divided by twice the variance. It is negative, indicating that fluctuations of adjacent terrace widths 
are anticorrelated. For Gruber-Mullins this covariance is ipso facto — 1. With the Grenoble formalism p3,nq] we find 



it to be -1/3 (NN) or -0.36. . . (all); with the Saclay formalism |21|, it is -0.33. In all these cases, the covariance is 
independent of A. In contrast, the exact covariance increases weakly in magnitude with A, from -0.31 at A— to 
-0.34 at A~ 2 Jnj, and presumably to -0.36 asymptotically. 
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Table 2: Compendium of experiments measuring the variances of terrace width distributions of vicinal systems. 



Vicinal 


T(K) 


a 2 


Q 


A 


Aw /Aq 


A w (eV-A) 


Pt(110)-(lx2) 


298 




2.2 


0.13 


— 


/3=? 


Cu (19,17,17) 


353 


0.122 


4.1 


2.2 


0.77 


0.005 


Si(lll) 


1173 


0.11 


3.8 


1.7 


0.96 


0.4 


Cu(l,l,13) 


348 


0.091 


4.8 


3.0 


1.27 


0.007 


Cu(ll,7,7) 


306 


0.085 


5.1 


4 


1.37 


0.004 


Cu(lll) 


313 


0.084 


5.0 


3.6 


1.39 


0.004 


Cu(lll) 


301 


0.073 


6.0 


6.0 


1.58 


0.006 


Ag(100) 


300 


0.073 


6.4 


6.9 


1.58 


/3=? 


Cu(l,l,19) 


320 


0.070 


6.7 


7.9 


1.64 


0.012 


Si(lll)-(7x7) 


1100 


0.068 


6.4 


7.0 


1.67 


0.7 


Si(lll)-(lxl)Br 


853 


0.068 


6.4 


7.0 


1.67 


0.1 


Si(lll)-Ga 


823 


0.068 


6.6 


7.6 


1.67 


1.8 


Si(lll)-Al V3 


1040 


0.058 


7.6 


10.5 


1.85 


2.2 


Cu(l,l,ll) 


300 


0.053 


8.7 


15 


1.95 


0.02 


Cu(l,l,13) 


285 


0.044 


10 


20 


2.12 


0.02 


Pt(lll) 


900 


0.020 


24 


135 


2.59 


6 


Si(113) rotated 


1200 


0.004 


124 3 


8xl0 3 


2.92 


(27±5)xl0 2 



Experimenters 



Swamy, Bertel |35[ 
Giesen |,|| 
Bermond, Metois Q 
Giesen [5,55] 
Giesen [5,53] 
Giesen [5,53] 
Giesen [5,53] 
P. Wang. . .Williams 
Giesen pipaj 
Williams |36[ 
X.-S. Wang, Willianisjp 

Fujita...Ichikawa |5q] 
Schwennicke. . .Williams (59) 
Barbier et a l. 2l"| 

Giesen ||5| 
Hahn. . .Kern |o| 
van Dijken,Zandvliet,Poelsema 



The estimate of A is obtained from the (normalized) variance using Eq. (^), except for the first-row entry, which 
is based on a direct fit using the 2-parameter CGWD. Aq is short for Aqm(NN)i the conventional estimate (cf. 
Table 7 of Ref. B). The column "Ay/ I Aq \ computed using Eq. jlO|), shows that for most systems the correction 
factor is of order half the ultimate asymptotic factor. For the copper entries, A\y is computed from Eq. (|l|), using 
P TSK = 2fc B T(a||/a^)sinh 2 (e/2fc B r) and kink energies e of 0.126 and 0.12 for vicinals to {100} [nominally (l,l,2n+l)] 
and {111}, respectively. In other cases, Ay/ is simply rescaled from Ref. S. In two cases, the values of the stiffness 
are not readily available. 
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